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Abstract 

We present evidence that a special class of gravitationally-coupled hidden sectors, 
in which conformal invariance is dynamically broken in a controlled way, exhibit the 
properties of dark energy. Such quantum field theories may appear while embedding the 
Standard Model in a more fundamental high energy theory. At late times, an effective 
dark energy field behaves similarly to an exponentially small cosmological constant while 
at early times its energy density partly tracks that of matter. 
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1 Introduction 



What is the structure of our Universe and what governs its evolution? Extensive cosmological 
observations have shown that the Universe is mostly cold, dark and accelerating (for a review 
see e.g. [I] and references therein). A large fraction of its current (and possibly future) energy 
density may be modeled by an effective negative-pressure fluid named " dark energy" . Despite 
intense investigation, the fundamental structure and origin of this fluid is, however, not well 
understood. 

A simple cosmological constant term in Einstein's equations provides perhaps the least com- 
plicated model of dark energy. The smallness of the observed cosmological constant poses a 
theoretical challenge however and no compelling argument for it has been formulated to date. 
Dynamical models for dark energy have also been proposed: quintessence, K-essence, fermion 
condensates, phantom, tachyonic, coupled dark matter and dark energy models, to name a few 
(see [TJ[2] for recent reviews). Quintessence models introduce scalar fields and postulate specific 
potentials that can track dark matter and lead to late time accelerated expansion. K-essence 
models, originally proposed in the context of inflation, modify the kinetic terms of scalar fields 
and also lead to accelerated expansion. For some time period the equation of state of the mod- 
ified scalar field is the one expected for dark energy. In another class of models, the (bosonic) 
long wavelength excitations of a fermion condensate play the role of a dark energy field. These 
excitations are self- interacting and at late times relax to an effective cosmological constant. In 
phantom field models, a scalar field exhibits a "wrong sign" kinetic term; even if the classical 
dynamics of such a field can be consistent, the fact that the energy of phantom fields is un- 
bounded below makes their vacuum unstable and their quantum theory not well defined unless 
some mechanism generates a lower bound on their energy. Coupled dark matter and dark 
energy models postulate certain interaction between (dark) matter and the dark energy field 
and attempt to explain both the observed acceleration as well as the observed similarity of the 
matter and dark energy densities today. Other possible explanations of late-time acceleration, 
not relying on additional (fundamental or effective) scalar fields have also been proposed - 
such as infrared modifications of general relativity. 

All these scenarios, proposed as a consequence of the experimental evidence for dark energy, 
can be made consistent with observations and constraints; their place in a more complete theory 
describing both the evolution of the universe as well as the known particle physics, remains 
to be clarified. The Standard Model of particle physics explains successfully all collider and 
other experimental data; it nevertheless has many theoretical shortcomings and it is presently 
seen only as a good effective theory that requires a high energy completion. Many possible 
extensions of the Standard Model have been proposed, addressing specific difficulties of the 
Standard Model. Some of these models, while motivated by particle physics issues, provide 
good candidates for dark matter. It seems natural to expect that a compelling solution to both 
the cosmological observations and particle physics issues should exist in the context of a more 
fundamental theory. 
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One of the proposed ways for embedding the Standard Model in string theory is by making 
it part of a quiver gauge theory - i.e. a theory with many gauge groups and fields transforming 
either in the adjoint representation or in the bi-fundamental representation of these groups. 
Such theories may be realized (either directly or holographically) [3] in terms of D-branes 
placed at various types of singular points. The additional gauge groups are either broken at a 
sufficiently high scale or are decoupled from the Standard Model by appropriate choices for the 
relevant coupling constants or other parameters concerning the details of the compactification 
to four dimensions. Such decoupled sectors nevertheless affect the observable physics through 
gravitational interactions — either through higher- dimension Plank-suppressed operators or by 
playing the role of dark matter (if the excitations of these sectors are sufficiently massive) or 
dark energy (if the excitations of these sectors are sufficiently light or even massless). It is 
therefore interesting to study such decoupled theories from this perspective. All scalar fields 
present in these theories transform in nontrivial representations of the gauge group making them 
not directly suitable for dark matter or dark energy candidates. In this paper we will discuss 
a class of quantum field theories which may appear in such a scenario. They are conformally 
invariant at the classical level but no longer so at the quantum level. Conformal invariance 
will be broken in an interesting way which will allow us to gain all-order information about the 
theory and identify the relevant gauge singlet (s) that can play the role of a dark energy field. 

Further motivation for analyzing theories of the type outlined above is the proposal jl] 
to embed the Standard Model itself in a conformal field theory. In this scenario, at some 
(high) energy scale particle physics is governed by a non-supersymmetric conformal field theory. 
Sometimes such theories have unexpected features (see e.g. [5], [B]) whose phenomenological 
consequences are worthwhile exploring. 

From a string theory perspective these theories may be realized in terms of a stack of D3 
branes probing non-supersymmetric A n _i orbifold singularities; they are quiver gauge theories 
with n nodes. Fields transforming in the adjoint representation are realized by strings connect- 
ing branes in the same stack while fields in the bifundamental representation are the lowest 
lying modes of strings connecting the images under the action of the orbifold group. One may 
extend such construction by adding further D7 branes; apart from matter in the fundamental 
representation of the quiver gauge group (which is realized by D3-D7 strings), such a con- 
struction naturally adds gravitationally-coupled gauge singlets. Regardless whether the quiver 
gauge theory is interpreted as a hidden sector (as we do in this paper) or if it is broken to the 
Standard Model, these additional fields may be suitable dark matter candidates. 

In the next section we will briefly describe a class of non-supersymmetric models in which 
conformal symmetry is broken dynamically; due to the special properties of their renormaliza- 
tion group flow, these theories can reach a regime in which they simplify dramatically, their 
entire dynamics being governed by a simple matrix scalar quantum field theory. In §3 we will 
describe a curved space version of this model and construct an equivalent model whose fields 
are singlets under the gauge symmetry of the original model and thus may play the role of 
a dark energy field. The cosmological implications of this model will be discussed in §4. We 
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present our conclusions and further comments in §5. 



2 Flat space models 

Quiver gauge theories have long been subject of active investigation. If the gauge groups are 
of SU(Ni) type, the action possesses a variety of double-trace terms required by the absence 
of the £7(1) factors. These double-trace terms renormalize independently of the single-trace 
part of the action. If the coupling constant of the single-trace part of the action is fixed, the 
double-trace coupling runs without bound at low energies implying that, in some regime, the 
double-trace operators dominate the dynamics of the theory. 

This scenario is realized in certain nonsupersymmetric theories in the multi-color limit, such 
as orbifold theories and non-supersymmetric /3-deformed theories j5], [6], [TJ. The former the- 
ories, which will be our main interest in the following, is obtained from 



C = -l(\r\N) Tr [ F^Ff" + D^D^j + A 2 [$ J , $ J ][$ / , $j] + xtfx + Xxi'l^i, x] 
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by an orbifold projection by a discrete group T C S77(4). Here the gauge group is SU(\T\N), 
the trace is normalized such that the identity element of the gauge group has trace equal to 

|r|jv, a 

~~ 9ym\^\N is the 't Hooft coupling constant, $ are six real scalar fields transforming 
in the vector representation of £0(6) and 'y 1 are the 5*0(6) Dirac matrices. As usual, covariant 
derivatives describe the coupling with the gauge field and are proportional to a single power of 
the coupling constant gyu- Four fermions, not shown in equation ([I]), transforming in the spinor 
representation of 5*0(6) and having Yukawa couplings with the scalars, complete the action. 
We will denote by g the representation of the elements of T in SU(\T\N), where they act by 
conjugation and by r g and R g the representation of T in the spinor and vector representation 
of 50(6), respectively. R g and r g are not unrelated; if one presents the vector representation 
of 50(6) as the 2-index antisymmetric tensor representation of 577(4), then R g = r g ®r g . The 
orbifold projection retains in the action ([I]) only the components of the original fields obeying 
the following relations: 

A^gA^ (f) 1 = R T g J g (f) J \' /^//v'V • (2) 

Moreover, the overall factor of T disappears and the 't Hooft coupling of the orbifolded theory 
is X = g YM N - 

Perhaps the simplest example of this construction is for r = Z 2 (i.e. \T\ = 2); its nontrivial 
element acts as 

g = diag(lAr, -Ijv) r g = -1 4 R g = 1 6 (3) 

inside the original SU (2N) gauge group, and on the spinor and vector representations of 50(6), 
respectively. The resulting theory [8] contains two SU(N) gauge groups, 6 scalar fields trans- 
forming in the adjoint representation of the first group, 6 scalar fields transforming in the adjoint 
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representation of the second group, 4 fermions transforming in the bifundamental (N, N) and 
4 fermions transforming in the bifundamental (N, N). 

The single-trace part of such orbifold theories has been analyzed in detail in [HI [ED] where 
it was shown that these terms are inherited from the parent theory; in the case of the action 
([1]) (together with its fermionic completion) they are finite. As we previously mentioned how- 
ever, the action contains additional double-trace terms [3 [6] which receive nontrivial infinite 
renormalization. For the simple T = Z 2 example they are 

6C = -f2oO IJ O u - f O 2 (4) 

with 

IJ = Tr[#$ 7 $ J ] - — Tr\g® K ® K ] O = Tr[0$ 7 $j] . (5) 

Such double-trace terms are generic in orbifold field theories and may always be traced to 
certain twisted auxiliary fields; they are not protected from renormalization by the general 
arguments of [HI [10] . As shown in [5] , they are in fact required by the renormalizability of the 
theory. The coupling constant / of a double-trace operator \0\ 2 runs at 1-loop as 

M^L=P f = vof 2 + 2 10 \f + aoX 2 , (6) 

where A = gy M N is the 't Hooft coupling, t> is the 2-point function coefficient of the operator 
O, 70 is the anomalous dimension of O and ao is the contribution of the single-trace Lagrangian 
to the /3-function of the double-trace coupling. If the discriminant 

D = 7^ - Av a > (7) 

then Pf has a nontrivial zero and / will flow to it. Some supersymmetry appears to be required 
[5] for such nontrivial zeroes. If supersymmetry is completely broken, as in the Z 2 model de- 
scribed above, there exists at least one double-trace coupling whose discriminant D is negative, 
leading to the following running coupling: 

f M ) = -12- + —tan — In — . (8) 

VO V \V V/Xr / / 

Here we defined b = \/—D and, for simplicity, the chosen boundary condition is f(fi r ) = 
—doX/vo- Thus, at weak 't Hooft coupling A, the double-trace parameter / varies very slowly 
for a wide range of scales. Still, it blows up towards positive infinity in the UV at M = fj, r e KVo ^ b ^' 
and reaches — oo in the IR at M = ii r e~ 1TV(> ^ bX \ z This unexpected runaway behavior of the 
coupling / was interpreted [5] in terms of a tachyonic instability of the string theory dual. This is 



3 One may expect that this singular behavior is softened by the 1/N corrections, which introduce a positive 
beta function for A, making it approach zero in the IR. We will restrict ourselves to the large N limit. 
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however not problematic as it has been argued from several different standpoints (TTJ H2] that 
the energy dynamically becomes bounded from below and the theory may flow to a nontrivial 
IR fixed point. 

The expectation following from this analysis is that along this flow, the twisted dimension- 
2 single-trace operators appearing in the double-trace operators that are generated quantum 
mechanically as illustrated in eq. (jlj) acquire nontrivial vacuum expectation values. To capture 
this effect and describe the consequences of this process it is therefore necessary to evaluate 
the effective potential for these operators; from the standpoint of the orbifold theory this is a 
multi-trace effective potential, the dimension-2n term having n traces. The explicit calculations 
of [5] show that the only effect of 1-loop corrections in the orbifold action without double-trace 
operators is to generate double-trace operators. All other possible terms are forbidden either 
by the analysis of [HJ [TO] or are derivative terms. The addition of double-trace operators to 
the tree-level action renders the theory renormalizable. It is not difficult to see that in this 
deformed theory, all n-trace terms generated at 1-loop level receive contributions solely from the 
double-trace deformation for all n < 3. Thus, for the purpose of finding the multi-trace effective 
potential, we may dispense with the large set of fields of the orbifold action and capture the 
dynamics of scalar fields by considering a much simpler action - that of a purely bosonic matrix 
scalar field theory with a specific quartic double-trace potential. To capture the effects of the 
ignored gluons and fermions all one needs to do is to renormalize the divergences and use as 
running coupling that of the orbifold theory (jSJ). It is important to notice that, as a twisted 
operator develops a vacuum expectation value, the potential of the undeformed orbifold theory 
becomes nonvanishing as well. This constant, which is positive and of the order of the square 
of the vacuum expectation value of the twisted operator, should also be added to the effective 
potential. 

While in general there will be several double-trace terms that have the features described 
above, one may focus on the one whose coupling constant runs fastest. This is the case because, 
as indicated by eq. ©, the runaway behavior - and therefore the appearance of nontrivial 
vacuum expectation value for a twisted operator - sets in first for this operator. Including the 
other double-trace terms may be treated as a perturbation and the running of their coupling 
constant will likely be affected by the emerging vacuum expectation value of the leading twisted 
operator. It is not always clear which orbifold group element generates the dominant double- 
trace term; we will generically denote this element by 7 G V C S77(|r|iV) as we will not need 
its detailed properties. It is however typically the case [6] that the operator whose coupling 
constant runs fastest is constructed out of SO (6) singlets. Thus, the action we will be interested 
in is just 

C = -NTr[rj^(d^)(d^)} - f(M) |Tr[ 7 $<l]| 2 (9) 

where without loss of generality we kept only two of the six scalar fields and have organized 
them into complex combinations. 

While simpler than the original orbifold action, the action (Q is still complicated mainly 
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due to the existence of the matrix degrees of freedom. Its vacuum structure and, as we will 
see in the next section, its gravitational effects, may be described in terms of a simpler theory 
containing a single scalar field to which it is completely equivalent. This is a standard approach: 
we introduce an auxiliary field that linearizes the quartic interaction 

L = -NrrTr[dfid v $] + jn^jW ~ P^h' 1 ®®] - ^Tr[ 7 $<5] (10) 

and construct the 1-particle-irreducible effective action of this field. This action contains the 
complete quantum information about the original theory regardless of the size of the expectation 
value of the auxiliary field (p (it is, in fact not clear a priori that such an expectation value is 
generated at ajlfl All higher-point potential interactions of the original scalar field $ that are 
generated at the quantum level are encoded in the higher-point interactions of the auxiliary 
field. This construction replaces the scalar field $ carrying gauge degrees of freedom by a single 
gauge-invariant effective field which may in principle be observable, albeit only gravitationally 
coupled with the Standard Model in the setup discussed here. 

The scalar theory (jUJ) encodes, through the scale dependence of the coupling constant f(M), 
the effects of the truncated fields on the dynamics of the scalar field appearing in the twisted 
operator. The reverse effects - i.e. the effects of the scalar field evolution on the gauge fields, 
fermions and the other scalars - are ignored. At sufficiently low energies, when the vacuum 
expectation value of the twisted operator (or, alternatively, for the auxiliary field p) becomes 
comparable to the energy scale, this back-reaction can no longer be ignored. Accounting for 
it is crucial for the arguments [7J [TTJ [12] that the theory flows to a nontrivial IR fixed point. 
Heuristically, around this energy scale, parts of the gauge fields, fermions and scalars become 
sufficiently massive and no longer contribute to the RG coefficients and a nontrivial beta func- 
tion for the gauge coupling is generated. Consequently, the running of the double-trace coupling 
06]) will be supplemented by the running of the gauge coupling and the resulting system is ex- 
pected to no longer have a runaway solution. On symmetry grounds, the fixed point theory 
is not expected to depend on the vacuum expectation value of the twisted operators, which 
is therefore irrelevant in a flat space setup. In curved space however, the vacuum expectation 
value sources gravitational field and its gravitational effects will not change qualitatively along 
the RG flow of the complete theory. 

In general, the orbifold field theory ([I])-© as well as the reduced model (j2]) are symmetric 
under the discrete transformation 7 — > where a; is a root of unity of the appropriate order to 
keep u'j an element of the orbifold group. Introduction of the auxiliary field does not break this 
symmetry which is preserved if the auxiliary field transforms as ip — > up, as may be readily seen 
from equation ffTU]) . This symmetry constrains the 1PI effective action to be solely a function 
of (pip; in particular, tadpole-like terms linear in <p or (p are forbidden. 

4 It should be mentioned that for the toy model in © the 1PI effective action is 1-loop exact. This is, 
however, not the case for the orbifold theory which is the motivation for this discussion and therefore we will 
not use this fact in the following. 
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The potential for the auxiliary field typically has a nontrivial minimum [7j , which translates 
into a nontrivial vacuum expectation value for the twisted operator Tr[7$$]. This can be 
translated into a nontrivial expectation value for the scalar fields. Such expectation values 
have two distinct effects: (a) they may yield a nonzero positive value for the tree-level potential 
and (b) at 1-loop they contribute a positive constant term to the potential. In flat space both 
contributions are, of course, irrelevant. The details of these contributions depend on the details 
of the orbifold group. 

In the next section we will place the model discussed here in curved space. One may find that 
a finite density of $ quanta is created (in a gauge- invariant configurations); to account for this 
possibility we also add a chemical potential. A systematic way of accounting for this is to add 
the number operator to the Hamiltonian which is then transformed to a Lagrangian framework. 
Since the auxiliary field ip linearized the quartic interaction, the behavior of the resulting field 
theory in the presence of a chemical potential is quite analogous to the zero temperature limit 
of an ideal relativistic Bose gas. Classically, the fields <3? are massless with the interaction terms 
playing the role of an effective mass. At the quantum level the auxiliary fields <p and (p acquire 
nontrivial vacuum expectation values, justifying this interpretation. 

Unlike systems of fermions, in the zero temperature limit the chemical potential does not enter 
explicitly the expression of the bosonic partition function. Rather, it is determined separately 
from the requirement that all particles are found in the ground state. If €q is the ground state 
energy (i.e. the effective mass of the $ particles) then the chemical potential is determined by 
the particle number density in the ground state: [f| 

1 kT 

n = ^fa-M) - 1 ^ = eo + — ^ e °- 

The standard relation between the grand potential and the partition function yields then the 
energy density as 

p = lim (—kT In Z) + an . (12) 

T-s-O 



5 In general it is not possible to give different interpretations to the modes of a field depending on their energy. 
However, modes with specified spatial momenta of any one field are a set of measure zero in the path integral. 
Due to the contribution of the integration measure, J dE J k 2 dk, modes with vanishing spatial momenta make 
vanishing contributions to the 1-loop grand potential in the continuum limit. In a theory placed in finite volume 
one may isolate the contribution of the lowest energy mode. Its contribution is volume independent and thus 
subleading for large/infinite spatial volume which is the case for an FRW universe or in Minkowski space. foThis 
is a standard approach which is used extensively in the treatment of Bose condensation; its advantage is that 
it allows particles to dynamically condense to the ground state or leave the condensate. However, the ground 
state condensate should not, in our case, be described as a classical field. From equation (|10p it is clear that 
the auxiliary field ip acts as effective mass for the fields $. As we will see in later sections, the auxiliary field 
acquires generically a nontrivial expectation value thus rendering $ effectively massive. This effective mass is 
the ground state energy. Any translational momenta for the ground state quanta will raise their energy and 
this take them out of the ground state. We should also note that, since they have vanishing momenta, in flat 
space the ground state particles do not contribute directly to the effective potential for the auxiliary field, their 
contribution being cancelled by the integration measure. 



8 



One can easily see that the first term is nothing but the zero temperature field theory partition 
function to be obtained from the action (Q. It is also the Too component of the auxiliary field 
stress tensor constructed from the quantum effective action. Thus, a simple way to include a 
finite $ density at zero temperature is to add a pressureless stress tensor with T 00 = fin to the 
stress tensor of the auxiliary field if. 

Clearly, such an addition manifestly breaks conformal invariance. From the standpoint of 
the orbifold theory, for which ([9]) is only a model capturing its essentia features, all violations 
of conformal invariance occur only quantum mechanically and are proportional to the beta- 
function of the double-trace operator (in a cosmological setting further effects proportional to 
the Hubble constant may also occur). For this reason one may expect that the particle- number 
density n is small; while we keep it throughout the general analysis, we will set it to zero when 
exploring the cosmological consequences of our model. 

In the next section we will discuss the curved space version of the quantum field theory 
detailed above and compute its partition function. It will turn out that, with a suitable choice 
of variables, the curved space calculation can be easily mapped to a flat space calculation. 



3 Curved space model 

The effective model we identified in the previous section captures the features of a quantum 
field theory whose only departure from scale invariance is encoded in the running of a certain 
double-trace coupling constant /. In coupling it to gravity we would like to preserve this 
feature; to this end, apart from the standard covariantization of index contraction, we will also 
add a conformal coupling, £ = 1/6, for the scalar fields $0 

L = y/=g~ (-Ng^Trid^dv®] -£i?Tr[$$] - /|Tr[ 7 $$]| 2 ) (13) 

where, as before, 7 is an element of the orbifold group T C SU(N) and R is the Ricci scalar. 

Similarly to the flat space theory, we proceed by linearizing the dependence on $ by intro- 
ducing auxiliary fields: 

L = (-Ng^Trld^dv®] - f + jCpip - ^Tr^" 1 ^] - ^Tr[ 7 $<5]^ . (14) 

Scale invariance implies that a rescaling of the metric g^ v can be absorbed by a field redefinition. 
In conformal time, an FRW metric is just that of Minkowski space up to a scale factor 01(77). 
The presence of the conformal coupling then implies that this scale factor can be eliminated by 
a field redefinition. Indeed, it is easy to see that in terms of the new fields 

<| = a(r/)$ = a(r]) 2 if (15) 

6 Such a coupling is required for the renormalizability of the theory in curved background; the value £ = 1/6 
ensures scale invariance. 
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the action is exactly that of the original Minkowsky space theory: 



1 T 



L = ~NTr[d fl ^d f1 ^} + jipip - ^Trpy- 1 ^*] - ^Tr[ 7 $$] 



(16) 



This is a reflection of the scale invariance of (1131) . 

The scalar fields $ appear quadratically. Therefore the 1PI effective potential may be easily 
computed by evaluating the determinant of the operator 



K B = - V ^d,d u + ^1~ 1 + ^1 



cff 



i In det Kb 



(17) 



On general grounds this determinant has both power-like and logarithmic UV divergences. 
These divergences may however be traced to the fact that we are focusing on the dynamics 
of the fields $ and are ignoring the other fields. If a quadratic divergence were present, it 
would generate a tadpole for (p\ in the complete orbifold theory such a term would suggest a 
perturbative generation of a mass term for some component of $ at two loops. This however 
is not expected to happen in the large N limit [9j [10] , consistent with the cancellation of terms 
linear in (p. We will therefore discard quadratic divergences from our effective model. 

With this clarifications, and if there exists a choice of fields tp and analytic continuation 
of momenta such that Kb does not have zero eigenvalues, the loop-induced potential for the 
auxiliary fields (p and cp is just 



cff 



(4tt) 



:Tr 



A 4 



with 



(19) 



Here we carried out the trace over momenta; the remaining trace is only over the gauge group 
indices. In the complete theory the logarithmic divergence present in equation (Tl8|) contributes 
to the running of the coupling constant / through the coefficient vo in the /3-function ([6]). After 
renormalization the cutoff scale A is replaced with the renormalization scale \x T . 
For a Z n orbifold, the Coleman- Weinberg potential (fl8|) is: 



v: 



N 2 



cff 



E 

fc=0 L 



In 4 



a k a k 



(20) 



In the calculation above (p maintains its role as a non-dynamical field as it was assumed 
to be constant throughout the calculation. It is however easy to see that derivative terms 
are also generated at the quantum level. They may be organized following the number of 
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fields on which the derivatives act. To this end one separates M. into a constant part and a 
position/momentum-dependent part and expands in the latter. The terms in which no deriva- 
tives act on the momentum-dependent part of Ai may be resummed and lead to a trivial shift 
of the constant part of M.. This is similar to doing perturbation theory around an arbitrary 
value of (p. While in general this would manifest an instability, this is not the case here because 
at tree level any value of ip is allowed (alternatively, because the potential for ip is generated at 
the same order as the derivative terms). With the same assumption as before, that momenta 
may be analytically continued such that for fixed fields Kb does not have zero eigenvalues, the 
first correction, containing two derivatives, is given byEI 

Tr[M(-q)M(q)} . (21) 

Certain care is necessary in the identification of the leading term in the momentum expansion, 
which depends on the details of the orbifold group and choice of group element 7. Regardless 
of its precise expression, the meaning of this term from the standpoint of the original theory is 
that of an effective contribution to the 4-point scalar amplitude. Following a strategy similar 
to the above it is not difficult to find higher-derivative corrections. We will however refrain 
from writing general expressions here. 



5L 



K 



-Tr 



p 2 + M 



M(-q) 



1 

(p + q) 2 + M 



M{q) 



- -Tr 



{p 2 + Mf 



3.1 An illustrative example 

The simplest example illustrating the discussion in the previous section is the Z2 orbifold theory, 
i.e. n = 2; mo = —ni\. The corresponding numerical coefficients are such that = 1 = a\. 
In this case the field (p may be chosen to be real. Repeating the discussion above and accounting 
for the fact that in this case M. always has a negative eigenvalue, we find that the effective 
potential is given by 

V ^ = + (4^( v ln iV^- 1 J • (22) 
The complete potential is therefore 

= J^o-J + ^ ■ (23) 

As discussed in the previous section between equations (jSJ) and (jUJ), we have added to the 
effective potential the effects of the scalar potential of the undeformed orbifold theory at the 
critical points tp of V e $. As discussed there, it is a constant of the order of the square of 
the vacuum expectation value of the double-trace operator - i.e. it must be proportional to 
(^o) 2 / 7 '• Co is a numerical coefficient of order unity, i 

7 The evaluation of this correction relies on the observation that [M(0), M.(q)] — 0. This holds as both M.(G) 
and M.{q) are constructed form mutually commuting matrices. 



11 



For this simple case the two- derivative term ( 12"T|) is just: 

N 



SL 



K 



6(4tt) 



1 &*(pd,i<p . 



(24) 



This will promote the auxiliary field <p to a dynamical field. Higher derivative terms have the 
structure 



(25) 



with numerical coefficients c n ; further terms, in which derivatives act on three or more fields, 
are also generated. We will neglect all such terms in the following, assuming that (p varies 
sufficiently slowly. 

To identify the action for the auxiliary field ip (and consequently its equation of state) it is 
necessary to restore the background FRW metric and also return to the comoving frame. It 
is therefore important to understand the behavior of the renormalization scale fi 2 . under this 
transformation. In curved space any cutoff or other scales should have a covariant interpreta- 
tion. Thus, fj% should be thought of as a square taken with the (inverse) metric; this implies 
that, as we restore the background FRW metric, we should also replace 



ji 2 — y a(rf) 2 ^ 



2. .2 

r 



The potential therefore becomes 



V Z2 = a{r)Y 



CO O 



4a(r]) 4 ip 2 



In 



- 1 



a(r/) 4 iV 2 /^ 

Thus, after identifying covariant quantities, the potential for (p is 

4<^ 2 



-9 V(<p) 



(26) 



(27) 



111 



- 1 



(2f 



Restoring the scale factor in the derivative terms also leads to the appearance of its derivatives 
(and thus of powers of curvature invariants) and of additional potential-like terms. They turn 
out to have simple covariant expressions, being expressible in terms of the Ricci scalar: 



5L 



N 



K 



6(4tt) 2 

N 
"6(4tt) 2 



g^dJa{rifv)dMri)\) 



a(r])*<p 
— "1 2 

-g -g* 1 df,(pd u (p + -<pR 
Lv 9 3 



(29) 
(30) 



The power of (p multiplying the Ricci scalar is determined by dimensional analysis; the ap- 
pearance of the last term is a direct consequence of the conformal invariance of the original 
model. 
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The derivative term may be brought to standard form by a simple field redefinition 



<P 



(4tt) 2 C 2 



AN 



(31) 



the new field ( has the canonical dimension of a scalar field. The complete (two-derivative) 
Lagrangian for this new field is 



SL 

v(0 



1 
2 



sTd^ + -C 2 R 



Co 2 9(4tt) 4 m 
f*° lGiV 2 / 



9(4tt) 2 



4iV 4 /i 4 



J eff 



2G 



R + 5L-V 



(32) 
(33) 
(34) 



where G is related to Newton's constant by G = 8irG^. The scale factor may also be restored 
in the higher-derivative terms (125j) . The Einstein-Hilbert term can be cast in canonical form 
by rescaling the metric by (1 — G( 2 /6). We will keep the action in its current form. 



4 Cosmological consequences 

In previous sections we argued that, if the Standard Model is embedded in string theory through 
a certain class of quiver gauge theories, then decoupled sectors of that theory yield effective 
actions of the type (134)) which interact only gravitationally with the usual matter fields. It is 
therefore interesting to explore the cosmological implications of such actions. 

While different in details, the model constructed above exhibits elements of classes of models 
discussed elsewhere. The existence of a field- dependent Newton's constant makes it similar to 
modified gravity modelsjf] As we will see in the following, the same effective Newton's constant 
(or the conformal coupling of the field Q will lead to the matter energy density acting as a 
source for £ and vice versso Finally, the existence of nontrivial derivative terms and potential 
terms for ( makes it similar to K-essence and quintessence models. While of course the field- 
dependent Newton's constant may be eliminated by a suitable Weyl rescaling of the metric, all 
the other features of the model survive this transformation. 

Let us proceed with analyzing the cosmological implications of the action ( |34|) ; as we will see, 
we will recover many desirable qualitative properties of cosmological parameters. We will also 
see that the effective cosmological constant will turn out to exhibit an exponential suppression 
compared to the renormalization scale. It is not difficult to see that the equation of motion for 

8 Let us note here that terms containing curvature invariants and matter fields are generic at higher order in 
perturbation theory. 

9 For similar reasons, such matter-dark energy couplings are also generic in modified gravity models. 
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the field ( and Einstein's equations are: 



6 



BV 


= 


(35) 




= G (^L + ^ + 67%, + T™) 


(36) 


T C 


= d,Cd v C-\g lx Af P daCd p C > +2V) 


(37) 


rpc 


= ^ 2 E, V + \ [g, u D(C 2 ) - V M V,C 2 ] 


(38) 


ST C 


= g^g^s/^x^ 


(39) 



where □ stands for the usual covariant Laplacian operator. In the equations above T*>, T c , 
57 < > and T m are the stress tensor of the £ field in the absence of the conformal coupling, the 
contribution of the conformal coupling to the stress tensor, the stress tensor due to a finite 
density of * quanta and the matter stress tensor, respectivelffl 

The second term in equation (|35l) is also a consequence of the conformal coupling. It is 
useful to separate as much as possible the evolution of £ from that of the scale factor a (or 
more generally the metric if one were interested in constructing the perturbation equations). 
Evaluating the Ricci scalar from the trace of Einstein's equations and replacing it in fl35l) leads 
to 

□C-^-(^l-^C 2 J^ + ^(^ + T-) = 0. (40) 

We notice here that both the matter density as well as a finite density of $ quanta source the 
evolution of (. The coupling between them, polynomial in (, is different from existing analysis 
of couplings between (dark) matter and dark energy, see e.g. [T31 dU ITS] . 

These equations simplify substantially for an isotropic and homogeneous universe with flat 
spatial slices, for which the metric takes the standard FRW forml^l 

ds 2 = -dt 2 + a{t) 2 dx 2 . (41) 



10 The equations (|36]) - ([39|) receive, in principle, contributions due to particle production due to the expansion 
of the universe. There is, however, no direct contribution to the equation of motion for the dark energy field 
C (I55|) and therefore the contribution of particle production on the evolution of C, is expected to be small. We 
will not include the effects of particle production. 

n In flat space, the vacuum expectation value of the field £ is the critical point of the effective potential. 
The curved space analog is the set of equations above in which one requires Q = to be constant. Unlike the 
flat space case discussed in footnote [5j these equations introduce a correlation between Co and the density uq 
of the ground state condensate of $ quanta. This is, however, a second order effect which arises due to the 
gravitational coupling and is therefore expected to be heavily suppressed. 

12 Absence of a solution with such an ansatz simply means that one or more of the assumptions — homogeneity, 
isotropy or flatness — should be relaxed. 
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The assumptions of isotropy and homogeneity also imply that ( = ((t) does not have any 
spatial dependence; then, the equations above simplify to: 

f + 3Hi + fcv + (l - f e) f + f C (-^ICI^ - r H" 3P») = (42) 

3H 2 = G(p c + p m ) (43) 

2- = -§( p C + 3p C + p m + 3p m) (44) 
a 3 

where p m = y — ^ G is the effective matter energy density and p m is the matter pressure, which 

vanishes for non-relativistic matter. We have also identified the energy density and pressure of 
the C fluid as: 

The equation of state for the ( fluid follows immediately: 



l + u; = - ■ ■ ^— . (46) 



3 



V + ±e + HCC+l2nJ^\(\i 



!(t) 2 



At very late times in an expanding universe, if £ asymptotes to a constant, it is easy to see 
that the equation above reduces to w — — 1. In general, we see however that the acceleration 
term ( can help bring the equation of state close to w — — 1 without necessarily requiring a 
small kinetic energy (or an extremely small mass for the ( field). This acceleration term is also 
a consequence of the conformal coupling in (|34"|) . 

It is curious to notice that the dependence of ST^ on the scale factor a — identifiable by the 
coefficient Uq — implies that a nonzero $ density acts from the standpoint of the Friedman 
equation like nontrivial negative spatial curvature. Unlike regular curvature contributions how- 
ever, this term does not affect the curvature of spatial slices and is therefore unconstrained by 
observations. 

In practice it is convenient to replace the acceleration equation with an equation describing 
the time evolution of the matter energy density. As usual, Einstein's equations fT43|) and ( f44|) 
encode the time evolution of the total energy density. The time evolution of the energy density 
of the C fluid is, of course, given by equations ([15]) and (J4"2|) . It thus follows that 



p m = -pC-3-(p c + p m +/+p m ) , (47) 



which we will use in place of (jHJ). If £ is fixed to a constant this equation has the usual solution 
p m oc a~ 3 . We will see that this behavior naturally occurs at very late times. 
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4.1 A special solution 



In general, the equations discussed above yield a nontrivial time-dependent profile for the field 
( which is sourced by the matter density and by the density of $ quanta populating the ground 
state. As discussed in §2, the existence of a finite density of $ quanta is a quantum mechanical 
effect (proportional to the /3-function of the double-trace operators) and is expected to be small. 
Thus, assuming no accidental enhancements, we will set no = 0. In the absence of matter and 
for uq = the evolution equations admit a simple solution with constant value of ( = (q. This 
solution is deformed nontrivially by a nonzero matter density, n / or even by a choice of 
initial conditions that set ((t = to) ^ Co or ((t = to) ^ 0. 
Setting ( = (o, the equations (142"]) . ( 115]) become 

2 4CoV(Co) + (l - ^oYZ = (48) 

C=Co 



3 " u ouy V 6*7 d( 

GV(Co 



3 (1 - f Co 2 ) 



H 2 (49) 



and, as discussed before, the matter density has the usual l/a(t) 3 behavior. The first term in 
equation (1481 is a consequence of the conformal coupling of (; in the absence of such a coupling 
Co is simply given by the position of the minimum of the potential. The second equation 
requires that its left hand side be positive. The solution to this equation is the standard 
scale factor in the presence of a cosmological constant and an effective Newton's constant 
G e ff = G/(l — G(q/6). As discussed before, the equation of state reduces to the usual w = — 1. 

Due to the nature of the equations ( 1421) . ( 143]) and (147]) it is natural to expect that, in the far 
future, any time- dependent profile for ( will asymptote to the solution described above. It is 
therefore interesting to discuss it in more detail. The solution of the counterpart of equation (1481 
in the absence of the conformal coupling may be readily obtained: 

Co = ±^=e^f. (50) 

Z7TV D 

We see that it is exponentially smaller than the renormalization scale /i r for small and negative 
values of the double-trace coupling /. Such values are allowed by the special properties of 
the renormalization group flow discussed in §2. As mentioned there, / runs very slowly if the 
't Hooft coupling is small; it is therefore natural to consider a double-trace coupling that is 
small (yet larger than the 't Hooft coupling) and fixed over a large range of scales. With this 
starting point, the solution to (148]) may be found as a series in G( : 



Co = Co l + ^4G n C • (51) 



n>l 



13 If the left-hand side of this equation were negative, the FRW ansatz for the metric no longer yields a solution 
to Einstein's equations; the appropriate solution becomes the anti-de-Sitter space. 
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Evaluating the left-hand side of equation ( I49p we find that the effective cosmological constant 
is, up to irrelevant numerical coefficients, 



Aeff 



(1 - f Co 2 ) 



V(Co) 



n>l 



) 



(52) 



i.e. it is exponentially small compared to the energy scale /i r that governs the dynamics of 
the theory. Clearly, the smallness of the effective cosmological constant is a consequence of the 
exponential factor which is independent of the constant Cq in equation (1231) . making the main 
qualitative features of our results independent of the precise value of c . 

4.2 More general solutions 

While simple to analyze, the solutions discussed above are not generic even in the absence 
of matter and of a finite density of $ quanta. A typical solution has a nontrivial time- 
dependence triggered by the initial conditions for the ( field which relaxes in the far future to 
the solution discussed in the previous subsection. Perhaps the simplest way to construct such 
solutions is to decouple the two equations f flSl) and (|43l by solving for the Hubble parameter 
from the latter and replacing it in the former. The positivity requirement (l4l?j) receives derivative 
corrections and becomes 



While decoupled, the resulting equations are lengthy and nonlinear; finding the complete time 
dependence of their solution requires a numerical approach. 

In the presence of a nontrivial matter density the three equations (j4"2l) . f|43|) and (j47j) cannot 
be decoupled; the system can however be solved numerically with sufficient accuracy. In Figure 
[I] we show such solutions for the field (, the scale factor a, the Hubble parameter H and the 
equation of state parameter w as a function of time, both in the absence of matter (left panel) 
and in the presence of a finite matter density (right panel). While we have not enforced a phys- 
ical normalization, these plots allow us to understand the qualitative behavior of cosmological 
evolution in the presence of the ( field, governed by the action (|34|) . 

We chose t = 1 to represent the present tinted and show it by the vertical line in all plots; 
as initial conditions for the scale factor we required that a(l) = 1. When matter density is 
included, we normalize it to the observed ratio p m (l)/p^(l) ~ 1/3. 

14 As in the previous section, we will continue to assume that uq = 0. 

15 The unit time interval is given by NG 1 ^ 2 — -/Vrpi an k, implying that for a realistic cosmology the number 
of colors N in the hidden sector has to be very large. This is consistent with the general discussion in earlier 
sections, where we kept only the leading N dependence in the orbifold theory. This also guarantees that 1/N 
effects cannot affect the running of the double-trace coupling, in particular the fact that it becomes negative 




Ga 2 - 



> . 



(53) 



(cf. eq. ©). 
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Figure 1: Numerical solutions for £, the scale factor a, the Hubble parameter H and the 
equation of state w. 
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Figure 2: The energy density of the ( fluid (red) and matter (blue) as a function of time. At 
early times the two curves are essentially parallel, showing that the ( energy density tracks the 
matter energy density Despite the oscillatory behavior of £ the energy densities are monotonic 
functions. 

As can be seen from the figure, the field ( oscillates at early times and eventually settles to 
a very small constant value. This behavior should indeed hold both in the presence and in the 
absence of matter, as the matter energy density diminishes at late times, implying that in both 
cases the solution should asymptote to the constant solution described in the previous section. 
Due to their coupling, a nontrivial matter density leads to an increase in early time amplitude 
of ( oscillations and in their frequency!^! At late times, the behavior of the system is that of 
a cosmological constant, with a constant Hubble parameter, exponentially growing scale factor 
and equation of state w = —1. At early times the equation of state of the ( fluid varies with 
time and w can reach positive values. In the presence of matter, the energy density in the ( 
fluid partially tracks the matter energy density. This is a direct consequence of the fact that 
the matter energy density sources (, as seen in equation ( H2"j) . As we previously discussed, this 
is due to the quantum mechanically generated conformal coupling 1 / 6R( 2 in ( 1341) . Fig. [2] shows 
the energy densities in the £ fluid (in red) and in matter (in blue). The position of the crossing 
point, shown in figure [2] at a{t c ) ~ (p m (l)/p ( '(l)) 1 / 3 , is determined by the absolute values of 
the matter and ( energy density. 

5 Conclusions 

In this paper we discussed a class of models that may appear if the Standard Model is embed- 
ded in a more fundamental theory through a quiver gauge theory. We assumed that, by some 

16 While not visible in the figure, it also increases the time necessary for £ to reach its asymptotic value. 




19 



mechanism, these models interact only gravitationally with regular matter and thus have only 
cosmological consequences. The flat space properties of these models were previously consid- 
ered: due to the special properties of their renormalization group flow, their dynamics is mainly 
governed by a certain double-trace operator. Making use of this observation we constructed a 
simpler effective matrix bosonic model that captures the essential features of the original. Its 
structure is such that it is completely equivalent to the quantum theory of a single effective 
scalar field. We computed the potential of this field which is 1-loop exact. Departing from 
the standard treatment of such models we also evaluated the quantum mechanically generated 
kinetic energy of this effective field and showed that it exhibits a conformal coupling with 
gravity. 

The example we discussed in detail suggests that, from a cosmological standpoint, such 
almost-conformally invariant "hidden sectors" mimic the expected properties of dark energy at 
late times, while partially tracking the matter energy density at early times. Asymptotically 
in the far future we found that the effective field becomes constant and thus its consequences 
are similar to those of a pure cosmological constant which is exponentially smaller than the 
renormalization scale of the theory. The exponential suppression is governed by the coupling 
constant of the original double-trace operator. The time evolution of the effective (or dark 
energy) field is in general nontrivial and sources an interesting behavior for the scale factor, 
Hubble parameter and equation of state. The value of the latter at present times is sufficiently 
close to 10 = —1. The conformal coupling plays an important role here, avoiding the usual 
extreme limits on the quintessence boson mass. In fact, the mass of the dark energy field is 
comparable to that of the far-future cosmological constant scale. 

While the model we discussed in detail predicts an exponentially small far-future cosmological 
constant, we cannot, of course, claim to have solved the cosmological constant problem. Indeed, 
the contribution to the cosmological constant of the zero-point energ y o f the Standard Model 
fields and perhaps other hidden sectors remains potentially uncanceledo Our discussion should 
be thought of pointing out an yet unexplored source of dark energy which by itself does not 
introduce any further problems at the quantum level. 

Since the hidden sector interacts only gravitationally with regular matter, it needs not be in 
thermal equilibrium with it. It is therefore possible to assume that its temperature vanishes, as 
we have done. Perturbative inclusion of a small temperature is not difficult: to leading order it 
corrects the effective potential though a temperature-dependent mass term for the dark energy 
field (. While the structure of the special solution discussed in §4 is unchanged, the effective 
cosmological constant receives a temperature dependent positive shift. It should be possible, 
though perhaps not straightforward, to analyze the effects of non-zero temperature beyond 
leading order. 

17 It is perhaps important to recall that various mechanisms for the cancellation of the standard model contri- 
bution to the cosmological constant in supersymmetry-breaking settings have been discussed in the literature - 
see e.g. [16] for a partial list of references. 
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Due to the special properties of the model, a finite density of hidden quanta contributes 
differently than regular matter to Einstein's equations: it is similar to a negative curvature 
term without actually changing the curvature of the spatial slices. It should be interesting to 
explore the consequences of such a condensate. 

While we formulated our analysis for the special case of quiver gauge theories which are 
conformal in the planar limit, the mechanism proposed here should hold in more general theories 
in the presence of double-trace operators. While we expect the general features to remain the 
same, the technical details will likely be different and model-dependent. In general however, 
due to absence of conformal invariance, there is no reason to introduce a conformal coupling 
for the scalar fields. It would be interesting to explore such more general models. 
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